function agg = kimball_agg(Y, eq, param, glob)

    agg = sum(eq.L.*upsilon(1/Y*eq.yf, param));

end